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The 'retarded' Green function for fields propagating on a Schwarzschild black hole spacetime 
possesses a branch cut on the complex frequency plane. Classically, the branch cut is important, 
for example, in order to fully determine the response of the black hole to a linear field perturbation. 
The branch cut is also useful for the calculation of the self-force on a point particle moving in the 
Schwarzschild background. In this paper we use techniques of analytic-continuation to the complex 
plane of the radial coordinate in order to calculate the branch cut contribution to the Green function 
in the limit of large imaginary frequency. It is expected that the contribution of this frequency regime 
to the perturbation response and to the self-force will be mostly for short time intervals. We also 
determine the highly-damped quasinormal mode frequencies for electromagnetic perturbations in 
Schwarzschild for the first time (previously only the leading imaginary part was known), which 
seem to have a 'deep connection' with the branch cut. We find that these frequencies behave like 

in i[£(e+l)]\ ■K^^''{l-i)[l{£+l)f , rj.u u- ^ ^ ■ , J 

ujin = — r ^-r — I ^ — , ; — hu (n I. 1 he highly-damped quasmormal modes 

are particularly interesting for theories of quantum gravity in that they are believed to probe the 
small scale structure of the spacetime. 



I. INTRODUCTION 



In the investigation of spin-field perturbations on a curved spacetime, the 'retarded' Green function plays a crucial 
role. For example, the full evolution in time of some initial data may be determined by integrating over space the 
'retarded' Green function convolved with the initial data. Also, the self-force acting on a point particle moving on 
a background spacetime may be calculated via an integration of the 'retarded' Green function over the whole past 
worldline of the particle [1] . Similarly in quantum field theory the covariant commutation relations are determined 
by the 'advanced' and 'retarded' Green functions. While in some contexts such as renormalization a knowledge of the 
short-distance behaviour of the Green function is sufficient, in other contexts the knowledge of the global behaviour 
of the Green function is important. 

The 'retarded' Green function for linear field perturbations of the Schwarzschild black hole spacetime can be 
expressed as a multipolc decomposition together with a Fourier integral over the frequencies just above the real axis. 
In the impressive work of 2 , Leaver deforms the frequency integral into a contour in the complex- frequency plane 
and he investigates three distinct contributions to the Green function coming from: (1) a high-frequency arc, (2) 
the poles (quasinormal modes) in the lower plane, and (3) a branch cut (BC) along the negative imaginary axis 
(NIA). Under a linear field perturbation, the high-frequency arc leads to a prompt response; the poles lead to the 
well-known quasinormal mode 'ringing' (this 'ringing' was observed for the first time in [3]); the BC leads to, at least, 
a power-law tail decay with time [4j[5]. While the quasinormal modes (QNMs) have been extensively investigated, 
most work on the BC has been limited to small frequency on the NIA (and usually also large radial coordinates), 
which is precisely the regime that gives the power-law tail - see, e.g. [HIS]- The BC, however, is also expected [5] to 
have a significant contribution at 'early' times, just after the start of the prompt response from the high-frequency 
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arc. In fact, at 'very early' times, i.e., before the prompt response, the separate contributions from the high-frequency 
arc, QNMs and BC might all separately diverge with their divergences cancelling each other out so that the complex- 
frequency contour integration is still valid. It is expected that these BC contributions at 'early' and 'very early' 
times come from the mid- and high-frequency regimes of the BC on the NIA, which might also bring about other 
yet-unsuspected contributions. To the best of our knowledge, these mid- and high-frequency regimes of the BC have 
only been investigated, respectively, in [3{9] and in Maassen van den Brink's jTOl, and it was done solely for the case of 
gravitational perturbations. In [10| . Maassen van den Brink investigated the BC in Schwarzschild for high-frequencies 
using a method based on an analytic continuation to the complex plane of the radial coordinate r. This method is 
mirrored on the method used in [HI [12] for the calculation of highly-damped QNMs, i.e., QNMs for large overtone 
index n, which therefore lie far down in the lower frequency plane. 

On the quantum side, highly-damped QNMs have attracted considerable attention since attempts [T51 - IT5] have been 
made at linking them to the area quantization of a black hole [121 [TT] . For example, in (TB] , the imaginary part of 
the highly-damped QNMs is shown to be related to the exponential redshift of the wave modes close to the horizon. 
In jl9j, highly-damped QNMs in Kerr spacetime have been interpreted as semiclassical bound states along a specific 
contour in the complex-r plane; they speculate that QNMs and another set of modes (the so-called totally-transmitted 
modes) correspond to different sets of microscopic degrees of freedom which, when they interact, produce Hawking 
radiation. Separately, in [20] they have shown that highly-damped QNMs probe the short length scale structure of 
the Schwarzschild black hole spacetime by calculating them in the context of a 'quantum-corrected' black hole. 

To set the current work in context we here give a succinct review of the -mostly- analytic results obtained in 
the literature for highly-damped QNM frequencies in the Schwarzschild spacetime for various spin-s fields. Motl 
and Neitzke [HI [T2| originally calculated analytically the leading-order of the highly-damped QNMs for s — Q and 
2 ( |12j also predicted the dependence of their next-to-leading order term on the 'angular momentum number' t). 
Subsequently, Maassen van den Brink [TO] for s = 2 and Musiri and Siopsis [H] for s = and 2 confirmed the 
leading-order behaviour and obtained the next-to-leading order highly-damped QNMs. In [21] [22] they also give the 
leading-order for the imaginary part -though not for the real part- of the QNM frequencies for s = 1/2, 1 and 3/2. 
Finally, in ^23j they obtain the leading-order and next-to-leading order QNMs for s = 1/2 and s — 5/2. To the best 
of our knowledge, the leading-order for the real part of the QNM frequencies for s = 1 and s = 3/2 remains unknown, 
and only the leading order of their imaginary parts is known. The case of electromagnetic (s = 1) perturbations 
is particularly interesting since the real part of the spin-1 highly-damped QNM frequencies have been expected to 
approach the NIA, thus hinting at an unexplored possible connection between the highly-damped QNMs and the BC. 
Indeed, in [51] [23] they show that the real part of the electromagnetic QNM frequencies can only go at most like 
for n — (X), and in 1241 they find numerical indications that it goes like n~^/^, with a coefhcient which is a 3rd order 
polynomial in i(£+ 1) (with undetermined polynomial coefficients). We will confirm this behaviour and determine the 



polynomial coefficients in Sec IV This behaviour is in contrast with the 0(1) behaviour of the corresponding QNM 
frequencies for s = and 2 and faster-decaying than the 0(n~^/^) for s = 1/2 and 5/2. See [2S| for a thorough and 
recent review of QNMs in different spacetimes; see also [26] . 

As mentioned, the contributions of the mid- and high-frequency regimes of the BC still remain largely unexplored. 
In the present paper we use the method in |10j (which was restricted to s = 2) to carry out an asymptotic calculation 
in the high-frequency regime of the BC contribution to the 'retarded' Green function in Schwarzschild for fields of spin 
s = (scalar) and 1 (electromagnetic). We also calculate the highly-damped QNM frequencies for electromagnetic 
perturbations for the first time in the literature (other than the leading-order of the imaginary part, which is known): 
we calculate these frequencies up to order n~^, that is, leading-order for the real part and up to two orders after 
leading-order for the imaginary part. The feature that these QNMs approach the NIA 'so fast' has made them 
very unyielding: like we show in this paper, in order to obtain the leading-order behaviour for large-frequency for 
spin-1 we are required to go up to two higher orders in perturbation theory than we are required for spins and 
2. Electromagnetic QNMs are increasingly important, since the detection of the electromagnetic counterpart of the 
gravitational wave emission by astrophysical sources might play a valuable role for localizing the source and obtaining 
further information about it |27| . For completeness, we also reproduce in this paper the results in [TU] for the case of 
s = 2 (gravitational) and in [3T] for the corresponding spin-0 QNM frequencies. 

One particular physical application of the BC contribution to the 'retarded' Green function is the calculation of 
the corresponding contribution to the response of the Schwarzschild black hole to some initial perturbation. In this 
paper we investigate the role played by the high-frequency part of both the BC and the QNM contributions to such a 
response in the case of a spin-s(=0,l,2) field sourced by: (1) initial data of compact support and (2) a non-compact 
Gaussian distribution in the 'tortoise' radial coordinate. 



In two other papers [5S1[5S], we will investigate, using completely different methods, the BC in Schwarzschild in 
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the mid-frequency and the small-frequency regimes along the NIA. In most of the figures in the present paper we 
plot the various quantities required for the calculation of the BC contribution to the Green function. In these figures 
we compare the high-frequency asymptotics obtained here with the results using the independent methods which we 
will present in |29j. The method in [291 has good convergence in a 'mid'-frequency regime which overlaps with the 
'high'-frequency regime, as shown here; it also overlaps with the 'small'-frequcncy regime, as we will show in [2S]. 

This paper is organized as follows: In Sec|ll]we give the general formulae for the BC contribution to the 'retarded' 
Green function. In Sec |III| we perform the asymptotic analysis of the BC Green function for high-frequency along 
the NIA. In Sec IV we calculate the highly-damped QNMs. In Sec|V]we investigate the contribution from the high- 
frequency regime of both the BC and the QNMs to the black hole response to an initial perturbation. In Sec VI we 
present some conclusions. 



II. BRANCH CUT CONTRIBUTION TO THE GREEN FUNCTION 



The 'retarded' Green function for linear field perturbations in the Schwarzschild spacetime can be expressed in 
terms of a multipole decomposition together with a Fourier transform as 



oo 

Gret{x,x') = ^(2£+ l)P,(cos7)Gr*(r,/; At), GY\r,r'-M) = — / 

Gi{r,r = , r> = max(r, r ), r< = mm(r,r ) 

W (uj) 

W{uj) = W[gi{r,uj),fi{r,uj)] ^ ge{r,u) — fe(r,uj) — 



00 + 2C 



did G£{r^ r'] uj)e 



— iujAt 



oo -fic 



At = t- t', 
(2.1) 



where c > 0, t and r are ~ respectively - the time and radial Schwarzschild coordinates, — r + rhlii{f — 1) is the 
so-called 'tortoise' radial coordinate, r/,, = 2M is the radius of the event horizon, M is the mass of the black hole 
and 7 is the angle between the spacetime points x and x' . Note that the physical region r G (r/j, oo) corresponds to 
r* G (— oo,(X)). A bar over a quantity indicates that the quantity has been made dimensionless via an appropriate 
factor of r/i, e.g., f = r/rh, w = uirh, t = t/rh, etc. The function W (uj) is the Wronskian of the two linearly 
independent solutions fi and gg of the following homogeneous radial ODE: 



dr 



Vir) ^ 1 - 



2 +uJ^-V{r) 



rh 



■)pi{r,uj) = 



(2.2) 



A = 



1) 



The parameter s denotes the spin of the field: s — 2 corresponds to axial ~ also called 'odd' - gravitational perturba- 
tions (in which case Eq.( |2.2[ ) becomes the Regge- Wheeler equation 30J), s = 1 to electromagnetic perturbations [5T| 
and s = to scalar perturbations [H [5]. The ODE Eq.(2.2| has two regular singular points at r = 0, rh and an 
irregular singular point at r = oo. These singularities in the ODE generally cause the radial solutions to have branch 
points at r = and in the complex-r plane. For real uj, the solutions fg and gi are uniquely determined by the 
boundary conditions: 



fi{r,uj) - e" 



— 7> 



-oo, 
-oo. 



(2.3) 



Strictly speaking, the limits in Eq.(2.3) should actually be r, 
them by 



=Foo, but following standard conventions we denote 
=Foo, with a certain abuse of language. We also have 



— 7> -1-00, 



(2.4) 



where A°^^ G C and G C are reflection and incidence coefficients, respectively. It is immediate that the Wronskian 



is equal to W (w) — —2iu!A™^. The boundary condition (2.3) also define fg and gg unambiguously for Im(aj) > when 
r, G K; these solutions are then defined for Im(a;) < by analytic continuation (see [9j for details on the region of 
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validity of these boundary conditions including r^, as well as uj). On the NIA, the boundary condition (2.3) becomes 
meaningless since it does not allow one to exclude the exponentially-decaying solution, e"'"*'^''* for and e"*"''* for gi. 
We are therefore motivated to analytically continue to the complex-r plane, and impose these boundary conditions in 
the regions Re(— zwr*) < for fi and Re(iajr*) < for gi, where they define the solutions uniquely (see also [5^I33| ). 

Leaver has shown that g^{r,u) has a branch cut (BC) which can be naturally taken to run along the line 
wr : — )■ — oo ■ i. If r > 0, then gi{r,oj) has a branch point at w = and a BC along the NIA, w : — — oo • i. 
In [33 they have shown that the existence of a BC in the complex-w plane is linked to the asymptotic behaviour 
of the radial potential: the exponentially-decaying Schwarzschild potential V{r) ~ e'"*~^ [A -I- 1 — s^] /r^ as — ^ — oo 
leads to poles in fg on the NIA (these poles, however, are cancelled out in Gg by the corresponding poles in W (w)), 
whereas its slower-than-exponential decay (with the exception of the centrifugal barrier) V{r) — X/r'^ ^ 2Xrh ln(f*)/r2 
as r* — )■ oo leads to the BC in gi. The function fg does not have a BC in the complex-w plane and so the BC 

of Gg{r,r';uj) along the NIA is due to the corresponding BC in ge. Note that the Wronskian W {lo), and therefore 
also A™^, inherits the BC on the NIA from ge (see, e.g., Eq.l28 [3l] and Eq.34 [2]). 

It is convenient to define i/ = iu £ C (and P = vvh), which is positive along the NIA. We also define AA{—iv) = 
Aj^{—iv) — A-{—iv) for any function A = A(lS) possessing a BC along the NIA, where A±{—iv) = lim£_j.o+ A{±t — iv), 
with V > Q. That is, AA is the discontinuity of A(u!) across the NIA. 

As mentioned above, fg has poles on the NIA; these lie at P = fc/2, Vfc e N (with the exception for s = 2 of 
the algebraically-special frequency v — vasi defined below, which is not a pole oi ft [71 [3]), and so it is convenient 
to define a new radial function: fg(r,—iv) = — sin {2TTD) ff {r, —iv). Accordingly, we define the Wronskian W (uj) = 
W ge{r,uj),fe(r,uj) . 



From the radial ODE Eq.(2.2) and the boundary conditions Eq.(2.3l, there follow the symmetries 



These symmetries lead to 



'ji(r,uj) = gg{r,-oj*) fi{r,uj) ^ fg{r,~uj*) if 
W{uj) = W* i-uj*), VweC 



ge-{r, -iv) = .g^+(r, -iv) if e M, 
W^{-iv) ^Wl{-iv), Vi^>0, 



(2.5) 



(2.6) 



so that the discontinuity of gi across the NIA is only in its imaginary part. Note also that although has a cut, 
\-^Tbj\ ^oes not have a cut. 
The BC contribution to the 'retarded' Green function is given by 

G^^(x,x') = ^(2^ + l)P^(cos7)Gf^(r,r';At), Gf ^(r, r'; Ai) = — / di. AG,(r, r'; -ii.)e-''^*, (2.7) 



The functions gt+ and gg- satisfy the same homogeneous, linear 2nd order differential equation (namely, Eq.(2.2)) 



Therefore, lS.gg will also satisfy this same differential equation and it may be expressed as a linear combination of 
the solutions gg(r,—iv) and gg{r, +ii>). In addition, both and gg^ satisfy the same boundary condition (2.3) as 
r — )> oo: gg±{r,—ii') ^ e"''*, and so Agg will not satisfy this boundary condition. To understand this, note that 



the behaviour for r — >■ oo of gg±(r^—iv) 



''* is dominant over that of gg{r,+iv) ~ e in Agg the boundary 



conditions determine that the dominant terms cancel so that Agg must be proportional to gg{r,+ii>). Finally, Agg 
must be purely-imaginary from Eq.(2.6) and, since we know that gg{r,+iiy) is real- valued, it must be [2j[7] that 



Agg{r, —iv) — iq{v)gg{r, +w), Vi^ > 



for some real- valued function q{v) which entirely characterises the BC 'strength'. 

fg{r,-iv)fg{r' ,-w) 



Using Eq.(2.8), we can express the BC modes as [31 [7] 

fg{r,~iv)fg{r' 



AGg{r,r'\-iv) = -2ivq{v) 



W+ {-iv)W^{-w) 



-2ii'q(y) - 



\W±{-iv) 



(2.8) 



(2.9) 



The second step is due to Eq.(2.6), and we include the double subindex in W± because A 



has no cut. 
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III. LARGE-// ASYMPTOTICS ON THE BRANCH CUT 

A. Method 

We essentially follow the method of Maassen van den Brink lOj, who calculates the high-frequency asymptotics of 
/f, gi+ and q{v) on the NIA for the case s = 2, to obtain the corresponding asymptotics for the cases s — {) and 1. 
For completeness, we will also include the results of [10] for s — 2. We now give a general description of the method. 

As described in the previous section it is convenient to impose the boundary conditions on and fg in the complex r 
or plane. The relationship between these two planes is complicated by the existence of a cut in the interrelationship 
r* = r + r/j ln(r/r?i — 1). Of special importance to our analysis are the anti-Stokes lines which for uj on the NIA are 
defined as the curves on the complex-r plane where Re(r'*) — 0. In this case there are four anti-Stokes lines emanating 
from near r = and they have slopes equal to ±1 (i.e., argr = ±7r/4, ±37r/4). The full structure is illustrated 
in Fig[l|a). We have two linearly-independent WKB asymptotic expansions for large-i>, ga{r^^iv) ~ e*'^''*; these 
expansions are valid for IfVi^l ^ 1 and away from the singularities r = and of the radial ODE. The importance 
of the anti-Stokes lines is that along them neither expansion, 5a('', +w) nor ga(^, dominates over the other. 

The contour chosen for g({r,—ii'), illustrated in Figjljb), is as follows. By analytic continuation, we can impose 
the boundary condition for ge{r, —iv) on the anti-Stokes line going to \r\ — >■ oo on the upper complex-r plane instead 
of imposing it for r — )■ oo [32, 33] • One can safely match gi{r, —iv) to ga{i', —iv) for \r\ — oo on that anti-Stokes line 
(i.e., at point A of FigQb)) since, there, (7a (^j —iv) does not dominate over ga(r, -\-iv). One can can then continue this 
expression from \r\ — >■ oo along the anti-Stokes line down to a region near r = with arg r = 37r/4. (i.e., down to point 
B of Figjljb)). However, one cannot continue this expression from point B to r > since one would have to cross 
into a region where ga{r^—iv) dominates over ga{r^-\-iv). Therefore, instead we introduce two linearly-independent 
solutions ipi, i = 1, 2, whose behaviour we can determine analytically and which we can match to the solutions we are 
interested in. Of course, we cannot determine these solutions exactly but we can obtain them in the limit of large D 
with fixed f\/l? and this is sufficient for our purpose. In this approximation, the solutions ipi are expressed in terms of 
special functions which one knows how to analytically continue from one anti-Stokes line to another. One then matches 
the solutions ipi to ga{r, —iv) ^ gi+{r, —iv) along argr — 37r/4 in a region of overlap, given by v"^/"^ ^ f ^ v^^/-^ 
(a region which includes point B of Figjljb)). Since one knows how to analytically continue ipi around r — with 
the matching done, one can analytically-continue gejf.{r, —iv) from the anti-Stokes line on argr = 37r/4 to another 
anti-Stokes line on argr = 7r/4 and have it expressed as a linear combination of ga{r,±iv). This linear combination 
is safely valid along the anti-Stokes line clockwise all the way to r* = (i.e., point C of Figjljb)). The coefficient of 
ga{r, —iv) in this linear combination is 1, and the coeflticient of ga{f, +iv) yields the asymptotics for the BC 'strength' 
q{v). From r* = 0, the linear combination is asymptotically valid on the physical line, either along r* € (— oo,0) 
or along r* e (0, oo), even though the coefficient of the subdominant solution in the corresponding region becomes 
meaningless. 

One can also find, using Frobenius method, two linearly-independent series expansions, ipi{r) and ip2{r), about 
r = 0. These expansions are valid for \ f'^v\ <C 1 and so there is no overlap with the region of validity of ga{r,iiv); 
the expansions ipf do however provide a check for the solutions ipi. 

The contour chosen for fi{r,—iv), illustrated in Figjljc), is as follows. First, from its boundary condition as 
r* — — oo we know that at r* = the solution fi(r^—iv) has to be a linear combination of ga{r,±iv) with the 
coefficient of ga{r,+iv) (which is the dominant WKB expansion on r^, S (— oo,0)) being equal to 1. This linear 
combination of ga{r,±iv) is valid anticlockwise along the anti-Stokes line from r* = up to argr = 7r/4. In order to 
obtain the coefficient of gair, —iv), one re-expresses this linear combination as a linear combination of ipi and ip2, on 
argr = 7r/4, using the matching described above. One knows how to analytically continue this expression onto the 
anti-Stokes line on argr = — 7r/4, and there it can be re-expressed as a new linear combination of ga{r,±iv). One 
can then continue this linear combination anticlockwise along the anti-Stokes line back to r* =0, thus yielding the 
asymptotic monodromy of fe around r = rh- The obtained asymptotic monodromy is then compared with the exact 
monodromy fg ((r — r/j)e^^^ —iv) = er'^'^^^ fi.{r — r^, —iv), which follows from the boundary condition on f( at the 
horizon. The comparison then yields the coefficient of ga{r, —iv) that we wanted. 

Let us note that in [32] they use a similar, but different, method to the one we use here. In [32] they make use of the 
Stokes phenomenon, instead of the solutions ipi, in order to continue from one anti-Stokes line to another. For that 
reason, they require the knowledge of the topology of the anti-Stokes lines very near r = 0, down to r = O (i?"^/^), 
where their WKB expansion breaks down - see Fig.l [52] • In our analysis, however, the WKB expansions never reach 
a region so near r — 0. See also [TT] for a description of similar contours followed. 
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Im(r) 




Rc(r) 



D 



FIG. 1: (a) In blue: schematic illustration of anti-Stokes lines (i.e., where Re(r,) — 0) on the complex-r plane in the case that 
cj is on the NIA in Schwarzschild. [See, e.g., Fig.l [31] for the topology of their anti-Stokes lines for large-frequency very near 
r = 0, where their method differs only slightly from the standard WKB analysis.] In dashed black lines: the region where 
Re(r*) < 0. (b) In red: the contour we follow in order to calculate the large-i? asymptotics for gi+. (c) In red: the contour 
followed in order to calculate fe. Point A corresponds to a point for |f| ^ 1 along the anti-Stokes line that stems from near 
r = at arg(r) = 37r/4; point B corresponds to a point on that same anti-Stokes line but lies within u^^^'^ <C <C i'^^^'^; point 
C corresponds to r, = 0; point D corresponds to a point with r > rh but r sufficiently far from rn so that the WKB expansions 
ga [r, -iziv) are valid. 



B. Radial solution gi{r,—ii/) 



Firstly, we find WKB asymptotic expansions for the solution of the radial ODE for ^ 1: 



51 W 



92ir) 



1 



dv 



vV{v) _ i [2\f +1- 



V{r) 1 
~A 8 



" , vV{v) 

dv 

V - I 



4^2 

n 2 



— 1) [Af + 1 - s2] - [2Af + 1 - s2] ' 



(3.1) 
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These expansions ga{r,uj) are valid away from the singularities r = (specifically, \T\fv\ ^ 1 is required) and Vh of 
the radial ODE and, they are such that neither expansion, Qair, +w) nor ^/a — ^^); dominates over the other along 
an anti-Stokes line. 

Secondly, we use the Frobenius method to find two linearly-independent solutions about r = of the radial ODE 
(re-expressed with r, not r,, as the independent variable). The characteristic exponents are: 1 ± s, and so they differ 
by an integer number when s = 0, 1, 2. We first find a power series solution tp[{r) about f = 0: 



s = 0: i^[{r) =f- Xf"^ + 0{f^) 

,F/ ^ 2 (A -2) 3 (A-2)(A-6) 4 ^,5, 



24 



s = 2 : ■^f (r) = f 



3 (A-6) . 



+ O(f^) 



(3.2) 



and the second, linearly independent solution V'f (^) is given by 

^ M 2 (2 + 2A-3A2) 
s = : (r) = [1 + 2A] f2 + 



+ 0{r'^) + ijj[ {r) In f 



s = l: VfW = l + Af + 0(f2)- yVfWlnr 



s = 2 : '02 = ^ + + + 0{r ) + 



2^ , 



(3.3) 



where 1/ = i^j^g = A(A — 2)/6 is the so-called algebraically special frequency (note that we use this term to refer both 
to Oas ctnd to loas = ~ivAs)- Note the exact monodromics: 



s = 0: 7/'f(re2") = Vf(r)-h27r#f(r) 
s = l: ^f(re2") =^/>f(r)-A27r#f(r) 



(3.4) 



s = 2: V2^(re2")=^f(r) + 



Vf(r) 



and (re^'^') = (r), Vs = 0, 1, 2. 

We now rewrite the ODE in terms of the independent variable t = f\fv and then group terms in different powers 
of v: 



(3.5) 



^-t|-(a^-l)-*^ L..2.3^^t^|-[A + .^-l], Z..-. 



Expanding for largc-;^, there arc two linearly-independent solutions '^j = '<l)f^ + V'l^\ * = 1)2, where '^j^'^-' denotes 



the two leading order solutions and '0-^^ the next-to-leading order ones: 

/■* 1 
Lo^f(t)=0, 4^\t) = du G^{t,u)Riiu), Ri{t)=^Ui,f'\t), Vz = l,2 



(3.6) 



G^{t,u) 



-v.(°)(i)vf (W)+V'f (i)4°^(«)" 









w 



'i,f\u),ijf\u)\=i,f\u) 



du 



This order will be sufficient for s = and 2. However, for s = 1 we will require the next order, ipi = ip^f^ +V'l^'' ) 



with: 



v,f)(i)=vr^(i)+vf)(t), i=i,2 



(3.7) 



^f"^(t) = / du G^{t,u)Rl"'\u), Rl"'\t) ^ ^L2^(t) 
Jo V 

)(i) = /* du G^{t,u)Rf\u\ Rf\t) ^ ^Ui>f\t) 
Jo vv 
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The leading order solutions can be found to be given by 



s = 0: <)(t) = _Jo(-), ^y^"it) 



t ^ ft 



ft 



(3.8) 



,(0).x e-"*/V^ , ft^\ 2sinh(tV2) .q), , e3W4^ 
s = l: = (^^ j = 4 ' W = ^*^i/2 (^^j = cosh (iV2) 



^/'t.. ft 



2i 



We have chosen the normalization constants so that, as indicated below, these solutions agree with the corresponding 
solutions tjjf. As for the quantities in the next-to-leading order solutions (and also the following order for s = 1), we 
find the following. For s = 0, 



W 



u 

V 



Ri{t) 



R2{t) 



(3.9) 



nt' 



\eY, (g)-(A-2i^)l^(^) 



4z/ 



For s = 1, 



W 



2u 



(3.10) 



Ri{t) 



V;«(i) 
R2{t) = 



2t [t^ cosh {ty2) - (A - 2t^) sinh {t'^/2)] 



y3/2 



4*3 cosh(iV2) + 3V^A e-*'/2erfi(*) - erf(t)e*'/2 

" 6^3/2 

t [i^ sinh {t^/2) - (A - 2t^) cosh {t^/2)] 



V 

' du 



1/2 



4°'(m)^2(u) V^A 



W 



2t 



12Asinh(tV2) + e-*'/2 {2 (e*'/^ - l) [t* - 3A] + 3X^t [e*'erf(i) + erfi(t)] } 



du ijf\u)R2{u) 



^ W 



12ti7V2 

For s = 1, we also required the following order. The first solution for s = 1 to following order is given by 
^(la) (^) ^ cosh(tV2) + (A - ^-^) sinh(tV2)] 

e-'-^ty/^e^ii + l)e*'erf(t) (A - 2t* - t^) - erfi(t) + 2^+1 {-2t^ + t^ + l)- 2t^ + t^l 



R'^'\t) 

e-Tt2 |a Ji2 ^ gt^ ^^2 + 3^ _ 3J _ ^2 J2t4 _ 13^2 ^ gt^' ^2^4 + I3f2 + 9) + 9] } 



3i>2 



^f«)(i) 



- cosh {t^/2) {-27eA 2A [chi(t2) - 21nt] } + 2Asinh (iV2) shi(t2) 



4i/2 



while the second sohition is given by 

-e-'^ty^i^^ii + l)e*'erf(t) {f+i- 2^^ - e) + crfi(t) + 2£2 + £ {~2t^ + t^ + l)- 2<4 + t^] | 



^ _ [(^ - cosh(tV2) ~ t2 sinh(iV2)] 



(3.12) 



4i> 

|a [-<2 ^ ^^2 + 3) + 3] _ ^2 \_2ti + 13^2 ^ (2<4 + I3i2 + 9) - o] } 



(0), 



(0)/ 




-Ini 
2 



W 



i^f\u)R'i'\u) A(A + 1) 



(A + 1)A 



hit 



(0), 



where erf(z) and erfi(z) — erf(iz)/i are, respectively, the error function and the imaginary error function, and chi(z) 
and shi(z) are, respectively, the hyperbolic cosine integral and the hyperbolic sine integral functions. Note that while 
shi(z) is an entire function of z € C, chi(z) has a branch cut along the negative real axis. Finally, the next-to-leading 
order solutions for 5 = 2 are given by 



W 



(0), 



Ri{u) -- 



TT 

4 



du 



Au 
v 



(3J3) 



^1 ^ 



•* du 

""2 
21 



^|2K-.„.^]n(^ 

/,(!) 



- Jl 

- lU Jq 

i?2(s) 4 

- iu'^Yn 



2i 
2i 



2i 



du 



R2{u) 



37r<3 j ' 



Note that the integrand for 11)2 given by Eq.(3.6) leads to an integral divergent at the lower limit, both for s = 1 



and s = 2. We therefore regularize it by integrating by parts and dropping the contribution at the lower limit which 
corresponds to adding an (infinite) multiple of the leading-order first solution, '0^''^ A similar divergence occurs for 
?/^2^"'* and i)'^^^ for s = 1, and we regularize them similarly. One can check that by taking an expansio n as t — > for 
tl'f\t) + V'i^''(i), * = 1, 2, both the leading order and the next-to-leading-order terms as r — 7> in Eqs.(3.2) and (3.3) 
are recovered for s = 0, 1 and 2. In the s 



ip''i\t) recovers the 0(r*) term of V^f ( r) in Eq.(3.2); similarly, the leading order term for t — > of 4'2^'{t) recovers 



1 case we have also checked that the leading order term for t 



of 



the 0(r'*lnr) term of ■02'(^) in Eq.(3.3) 
Using the formulae from |37j , 

Yk (ze™'^^ 

valid Vz, /i G 



(2)/ 



) = e"A<" (z) 

) = 6-"'''''% (z) + 2i sin(m/x7r) cot(/i7r) (z) 



(3.14) 



(-1)™^- [Yk (z) + 2imJk (z)] , 
and Vm, A; G Z, we obtain the following analytic continuations for the various terms in ipi. For s = 0: 

(3.15) 



<)(*i)=zVr(i), 



V.r(^^) = ^4"^(^)-|<^(^), 



?4^\t) 
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For s = 1: 



































TTZA (0) 












TTiX(X + 

^ 4 

















Note the terms with 4i'\t) in ^^2^°^' {it)^ i^2^'''{it) ^-nd tp'^' [it] arising from the regularization terms. 
For s = 2: 



/,(2a). 



,,(2b)/ 



<)(z<) = -z<)(t), 



(3.16) 



(3.17) 



It is also easy to check that, along argi = 7r/4, we have 



5 = 2: Vi"^e"'"^/' e 



(3.18) 



Next, we carry out an asymptotic series expansion for large-|t| along argf — 7r/4 of il^f^ and 4i^^ (and V'P^ for s = 1). 
In the corresponding expressions which we write next, we only show the t'^ terms in each of these large-|t| expansions. 
Using Eqs.(3.18l, we can then write for s = 0: 



JV2 



ai 



^/e*V2 



4y/D { 



-i + 



Pi 



+ e" 



(3.19) 



for s = 1: 





^ 0L2 




iyl/2 


V . 


+ ( 


fix 




-A 


iyl/2 


r 

V 



+ e-*V2 


-1 - 


PI' 


1 


















* ' 




V2 


-1 - 


^ + 




} 













argt — 7r/4 



(3.20) 



-tV2 



arg i = TT /4 



and for s — 2: 



1^ 



|e*V2 


H 


"1 ' 





















(3.21) 





-i 4 


A' 




1 + 


ft*" 















+ 0(t-2)+c> +o(/>-i)l argt = V4 



which is valid for \ <^t <^ v^l^ We can obtain the corresponding asymptotics f or tpi, 2{t) at argi = 37r/4 by using 
Eqs.(3.19)-(3.21 ) along argi = 7r/4 together with the analytic continuation Eqs.( 3.15| -(3.17). The dominant term 

of such asymptotics at argt 37r/4 (i.e., the term with e* must agree with the corresponding term in '>pi^2{it), 
obtained from Eqs.(3.19)-(3.21 ) for argt = tt/A. Similarly, we can take a linear combination of the asymptotics for 
"01.2 (i) at argt = 37r/4 such that the dominant term is zero, and then require that the remaining, subdominant term 
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(i.e., the term with e * agrees with the corresponding term in the hnear combination of ipi{it) and ip2{it), which 
are obtained from Eqs.(3.19)-(3.21 ) for argi = ir/A. The result of these requirements is: 



s = 0&2: aiGlR, /3i = -(2 + i)ai 

s = 1 : ai, a2 e M, /3i = -ai, Im(^2) = "i- 



(3.22) 



The actual values of ai can be obtained by directly expanding ip[^\t) in Eq.(3.6l for large \t\ with aigt — 7r/4. We 
obtain 



s — : ai — — 
s = 1 : ai = — 



487r3/2 



[1 + 3A1 



(3.23) 



From the boundary condition Eq.(2.3) and the leading-order behaviour of Eq.(3.1) it follows that 

gt+{r,-ii') ^ ga{r,-iv), i> > 1, for | argt - 7r| < 37r/4, 



(3.24) 



where we are only neglecting exponentially-small corrections. Performing now a power series in t/v^ with & > 1/6 
for ga{r,±iv)e^™'^e^*' 1"^ (note that we must replace er^^'"^ — e^*'^'^ when arg(f — 1) e (0,7r] arg(7~ — 1) e (— tt, 0), 
since this exponential comes from the ln(f — 1) in = r^(r)), we obtain 



1 + O^Tdi^ 
1 + O^Tdi- 



if arg(f - 1) e (0, tt) and t < P^/^ (3.25) 
if arg(f - 1) e (-TT, 0) and t < i^^/^ 



We will only require the coefficient d\ for the case s = 1, and in that case it is d\ = — A/12. We can then express the 
right-hand side of Eq.(3.25) as a linear combination of "01 (i) and "02(0 fo'' ^rgt = 37r/4, by using the asymptotics of 
Eqs.(3.19)-(3.21 ) at argt = 7r/4 together with the analytic continuation Eqs.( |3.15| -(3.17). The result is 



s = 



3»+ ^''j!^"^ +0(z>-^) 



^i(i)- 



2v^ 



-1 



ai 



Mt) 



P _ ^ A + 12mf-6(/32+/?2*) + 36A ^ ^ f j^-3/2\ 
2 yjl) 12D V / 

-A-12mf + 12a2 + 6(/32-/32*) ^ /-3/2\ 
12i? V ) 



s = 2 



+ {r, -zz/)e-"^ - ^ f3zi^3/2 ^ (2 - 3i)aiD + O (d^^^)] 0i(t) + ^ 

O L V /J -v/TT 



1 ai 



O 



which are valid Vargi G [0,27r), by analytic continuation. Usi ng th e asymptotics of Eqs.(3.19)-(3.21 ) for V'1,2 and 
rewriting the factors e^* in terms of the ga{r,±ii') using Eq.(3.25), we find 



s = 0, 2 : g,+ ir, ^iv) ^ g,{r, -iv) + (_l)i+^/22e2-- 
s^l: ge+{r, -w) ^[l + O (i^^^)] ga{r, -iv) + 

9a{r, +ii^) 



ai 



ga{r,+w) 



(3.26) 



V V \ ) 



The WKB expansion (7a(r, -|-ii/) is dominant over gaij, ^iv) in the region bounded by anti-Stokes lines which contains 
r = rfi, and so its coefficient can be trusted there. 

In Figj2]we plot the high-frequency asymptotics for of Eq .(3.26l together with Eq.(3.1 1 (we only include a token 
plot, for s = 0, since g^^ does not appear explicitly in Eq.(2.9l). The figure shows that these asymptotics agree with 
the completely independent calculation of . 
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FIG. 2: Log-plots of as a function of vM for s = 0, £ — 1, r — 2.8M. In blue: solution found using the method of [23]. In 
green: large-i? asymptotics of Eq. (3.261. (a) |Re((;f+)j, (b) ilm(5r£+)|, (c) 'zoomed-in' version of (b). 



C. BC 'strength' q{iy) 



It follows from Eqs.(2.8) and (3.26) that Agi{r, —iv) = 2ilm gi{r, ~iv) = iq{v)gi{r, +iv) with the following large-i/ 
asymptotics for the BC 'strength': 



s = 0,2 : 
s = 1 : 



Oil 



cos{2iTv) H — ■j= sin(27ri') 



0\ 



q{v) ^Aai 



sin(27ri/) ai cos(27r^) 



O 



Alternatively, considering q{v) in modulus- argument form we can rewrite Eq.(3.27l as 

ai 



s = 0,2 
s = 1 : 



q{v) ^ (-1)1+^/24 
1 



cos 27ri/ — 



q{v) ^ Aai 



sin l-KV 



since gairjiiv) G M and ga{r,+ii>) gi{r, neglecting exponentially-small terms. 



(3.27) 



(3.28) 



Figs IMS^ show that the large-i/ asymptotics of Eq.(3.27l match with a calculation of q{h') using the independent 



method of [29, (where we will show more clearly that the curve corresponding to the latter method agrees with Figs. 2 
in [7j and [5]). The asymptotics of Eq.(3.28l seem to do slightly better than those of Eq.(3.27) for s — 0, slightly 
worse for s = 2 while, for s = 1, they seem to do better in the phase but worse in the amplitude. 
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FIG. 3: BC 'strength' of Eq.( |2.8| as a function of v for s = 0, £ = 1. (a) Green curve: large-i> asymptotics of Eq.( |3.28l l; 
blue curve: q(y) calculated using the method of '29! via Eq.(2.8l with the choice of value r = 5M. (b) Same as (a), where now 
we also include in orange the asymptotics of Eq.(3.27|. 




FIG. 4: q(y'] as a function of v for s = 1, ^ = 1. See caption in Figjsjfor a description of the curves. 



D. Radial solution ji(r,—iv) 



From the boundary condition Eq.(2.3) and the relationship = r*(r), the following exact nionodromy around the 
horizon follows: 



(3.29) 



For r near (but far away enough from it so that the WKB expansions ga are valid), it is /^(r, —iv^ ^ e ^"^^ ^ ga{r, iv) 
to leading order as P ^ 1 and, since Qai'r, iv) is there dominant over Qai^, — w)^ we must have 



/^(r, -iv) - ga{r, iv) + c{v)ga{r, -iv) 



(3.30) 



for some function c(j^), where in the coefficient of gairyiv) we are potentially neglecting exponentially-small terms. 
Eq.(3.30) can be continued from the region of r 'near' to the point = and then anticlockwise along the 



anti-Stokes line up to arg(r) — 7r/4. From Eq.(3.25), we then have that 



t{r,-iv) ^ e e ' 



di 



c\y)e e 



di 



(3.31) 
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FIG. 5: q^u) as a function of u for s = 2, £ — 2. See caption in Fig|3] for a description of the curves. Cf. Fig. 2 [7] (also 
Fig.2 ^). 



for l/^/D < |f| < 1 with arg(r) = 7r/4. From Eqs.( |3.19| -( |3.2l| we find that 



0: fi{r,-iv)^ 

^ r 



VP 



s = 1 : /f(r, -ii/) 

^ ' —ITTU 



fi{r,-w) 
2 



_^ iai ia\ — 'Re {(32) + A 
\fl> V 

1 + + - "2 - ^Im (^2) ' 

, (i-2)ai 



^-1 + 



1 



ai iaf + Re (/32) - ^ 
\J1> V 



](- 




4) 







(3.32) 



-c(i^)e' 



i + 



(i + 2)ai 



^^1 + 





"1 ' 


+ c(z/)e"'' 


'1 + ^' 


}V'2 













This expression has been obtained with all functions /f, t/'i, '!/'2 evaluated for arg(r) = 7r /4, bu t it is valid Varg(r). In 
particular, we can express V'l and V'2 on arg(r) — — 7r/4 in terms of e^* by using Eqs.(3.19 1-(3.21 1 on arg(r) = 7r/4 
together with the analytic continuation Eqs.(3.15 1-( 3.17"| . We can match this expression to a linear combination of 
(7a(±ii/) via Eq.(3.1). This linear combination can then be continued anticlockwise along the anti-Stokes line all the 



way back to = 0, thus yielding the asymptotic monodromy: 

s = : /, ((r - rh)e^-\-iv) = e-2--gJr ~ iv) + 
s = l: j,[{r-Tn)e^^\-iv) = 



-2i 1 - 



gair~rh,~iv), 



(3.33) 



0\v 



2 : 



fi{{r-rh)e'^\-iu) 



-3/2 



c{v)e 



l + 0[v 



-3/2 



ga{r - Th.iv) + 



2i 1 



gair - Th, -iv), 



ga{r - rh,-iu), 
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Finally, comparing with the exact monodromy Eq.(3.29) and using Eq.(3.30), we obtain 
s = 0,2 : 



:(.) ^ (-1)-/^ : 2'^^^ . 



(3.34) 



1 : 



2iai 2ial 



O 



as well as, for 5 = 1 and from the coefficient of ga{r — r/^, ii^)^ 



-27rii/ 



l + 0(i>-^/2) +c(,y)0(i?-^/") 



-2i sin(27rz/) 



3/2\ 



(3.35) 



Note that the terms 0{D~^^^^) in Eq.(3.35l may actually decrease faster than v^^l"^ and they do not necessarily have 
to decrease both at the same rate. In particular, this implies that, e.g., c(y) may increase/decrease with T> as long as 
the 0{v~^l^) multiplying it makes up for it by decreasing/increasing faster than the other 0(z/~^/^). 

In Fig[6]we p lot = — sin [2 'kv) fi{r, —iv) as a function of Mv: it shows that the large-P asymptotics of Eq.(3.30) 



(together with (3.34) and (3.1 1) overlap with a calculation of fi using the so-called Jaffe series p9l 134 




FIG. 6: Log-plot of l/f] as a function of vM. The green and orange curves are the large-z/ asymptotics of Eq.( 3.30| (together 
with (3.341 and (|3.1[)) for r* — Q.2M and r, — 0.4M respectively. The blue and red curves are respectively calculated using 
the Jaffe series |29l I34| for r, = 0.2M and r. = QAM. (a) For s = 0, £ = 1. (b) For s = 1, £ = 1. (c) For s = 2, ^ = 2; note 
that both curves for the Jaffe series give a zero value for fi aX u = i>AS because of the definition of fi and the fact that, for 
s — 2, vas is not a pole. 
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E. Wronskian 



Wronskian on the 4th quadrant of the complex-w plane infinitesimaUy close to the NIA 



Using Eq.(3.26) for gg+ and Eqs.(3.30), (3.34) and (3.11 for fg, we find the following large-j/ asymptotics for the 

(3.36) 



s = 0,2 
s = 1 



W[ga{-w),ga{w)] 



sin(27ri/) 



sin {2tiv) + 2ie 



1 - 



(1 + i)ai 



21. A™ + = W+ {-iv) ^ W[ga{-iv).ga{ii^)] U + O (v-^'^) + 



2ai 2ial 



+ 0(9-'/') 



To leading order for s — Q and s = 2 it yields W^{—iv) ^ 2i. [1 — 2icot{2'Kv)]. This agrees with Eq.2.17 in (12] : 



lim 



>o+ V e 



1 



-AlTUJ 



[I + 2cos(7rs)] 



I. > 1, 



(3.37) 



u) — — e—i'i 



on the 3rd quadrant infinitesimaUy close to the NIA, after using the symmetry Eq.(2.6) 



Similarly to the 'trick' used to go from Eq.(3.27) to Eq.(3.28l, we can rewrite Eq.(3.36) as 



s = 0, 2 : W+ {-iv) 



W[ga{-iv),ga{iv)] 



sin {2'Kv) + 2ie^'"^^ 



sm{2'Ku) 

= 1: W+{-iv)^W[ga{~ii^),ga{ii^)]\l + 0{u-'"^] - 



(3.38) 



2ai 2x»p-<l+^ 
-e ^ 



V sin {2'Kv) 



where we have used Eq.( 3.34[ ) in the s = 1 case 



In Figs 



we plot w) = W gi^{r,—ii'),fe{r,—ii') as a function of Mv and we show that the large-z. 



asymptotics of Eq.(3.36) agree with the calculation using the method in 



F. Green function modes 



We can finally give an asymptotic expression for l arge- f or the BC discon tinuity of the 'retarded' Green function 
modes. From Eq.(2.9l and the leading orders of Eqs.(3.36), (3.30) and (3.27), we find 



AGi{r, r']~w) 
AGi{r, r']-w) 



(-l)"/22i cos(27rP) 
V [1 + 3cos2(27rj>)] 

-y/TTiI{i + 1) sui{2'Kv)rh 
P372 



(-l)''/2e'^''- +sin(27ri?)e~'"^*l U-iy'^e"'''' + sin(27rP) 



V^ljt + 1) 
2i^i/2 sin(27ri>) ' 



+ 1) 
2i?i/2 sin(27ri?)' 



, s = 0,2 
(3.39) 

s = 1 

- The 

divergence in the BC when At < |r*| + |r^| was to be expected, since the QNM series also seems to diverge at these 
'very early' times i2ji6j, and one would expect the divergences in the different contributions to the Green function to 
cancel each other out (see, e.g., [35] for the case of a radial potential which does not lead to a BC, and where the 
divergences of the QNM and high-frequency arc contributions at these 'very early' times cancel out). In Figs 10 - 11 
we plot the lavge-D asymptotics of AGe{r,r' ;uj) together with these modes calculated with the method in [29] . 



This shows the convergence of the {/-integral in the upper limit of integration in Eq.(2.7) when At > \r, 
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FIG. 7: Real part, imaginary part and absolute value of W+ as functions of fM. The asymptotic expressions in Eq.(3.36l 
are plotted as dashed black, orange and green curves, corresponding to the real part, imaginary part and absolute value of 
W+, respectively. The calculation using the method in [22] (for this, the value r = 2.8 Ad has been used to calculate the radial 
solutions gi+ and fi) is plotted as continuous blue, red and brown curves, corresponding to the real part, imaginar y par t and 
absolute value of W+, respectivel y. Th is plot is for the values s = 0, £ = 1. (a) Asymptotic expressions using Eq.(3.36l. (b) 
Asymptotic expressions using Eq. 13.381. The odd non-oscillatory intervals in the dashed green curves are just a visual artifact 
of the computational software program that we used for plotting the curves. 
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FIG. 8: Real part, imaginary part and absolute value of W+ as functions of vM. This plot is for the values s = 1, £ = 1. See 
caption in Fig[7|a) for a description of the curves. Note that, in this case, the method in [2^ (corresponding to the blue curve) 
cannot reach higher values of vAd. 



150 




t -'"•I 

FIG. 9: Real part, imaginary part and absolute value of W+ as functions of vM. This plot is for the values s = 2, £ = 2. See 
caption in Fig[7|a) for a description of the curves. Note that the blue and black curves are mostly overlapping, and so are the 
orange and red curves. The curves obtained with the method in [29] have a zero &t v = vas due to a corresponding zero in ft. 
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FIG. 10: BC Green function modes AGf(r(r« = 0.4), r'(r, = 0.2);— iu) of Eq.(2.9l as a function of v. The green curves 
correspond to the large-;? asymptotics of AGe using Eqs.(3.27l, (3.361 and (3.301 (so not just the leading-order Eq.(3^^[}. The 
blue curves have been obtained with the method in [55] (for which we have used fi instead of fe everywhere in Eq.(2.9l and 



we have chosen to calculate the radial functions at r = 5M and at r = 2.8AI for the calculation of, respectively, q{i') and the 
Wronskian). (a) For s = and £ — 1. (b) For s = 1 and 1=1. (c) For s = 2 and £ — 2 (the 'particular' behaviour around 
i/ = 4 is due to this value being that of the algebraically-special frequency). 



IV. HIGHLY-DAMPED QNMS 



QNMs are poles of the 'retarded' Green function in the lower frequency plane. By requiring that W+ = —IvA^l^^ 



in Eq.(3.36) be zero, we obtain the highly-damped QNM frequencies 

/2r4(l/4),_, .,3A + 1 ^/l 



In 3 ( n 1 

. = 0: ^,„^_- - + 



1447r5/2 



-(1 + 0- 



n 



ln3 fn l\ y2r4(l/4),, ,A-1 ^/l 



(4.1) 



where n is the so-called overtone index. These expressions agree with |21j (for s = 2 we are just including the result 
in [lOj for completeness). 

For s = 1, [211 [23] show that the 0(1) and 0(n^^/^) terms in the QNM frequencies are zero as n — J' oo. In [24] 
they find numerical indications that for s — 1 the highly-damped quasinormal modes go like 



in asA^ -I- a2A^ + fliA 



n — )• oo 



(4.2) 



with undetermined polynomial coefficients ai, a2 and 03. For s = 1, from Eq.(3.36) we would find the QNM condition 
to be 



1 + O + 



"2ai 2ial 



O 



c(P)e2'^*'' = 



(4.3) 
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FIG. 11: Similar to Figs |10| but using Eg. ( [3.38] ) instead of Eq.( |3.36l l. It shows a slight improvement with respect to Figs|10[ 

Let us try first with t he asymptotic expression 27rP ^ mr + a jr? ^"^ , for some undetermined coefficient a. We then 
obtain, from Eq.(3.34), that c(y) — 0{n) — 0{P), but in Eq.(4.3) we obtain a leading order 0{y/ri), which cannot 
be canc elled. Trying then with 27rP ~ nn + a + . . . or 27ri/ ~ nir + a^/n + . . . yields c(i^) = 0(n~^/^) but the 1 in 
Eq.(4.3) cannot be cancelled. Trying next with ^-kv ~ n7r + a/-y/n + . . . yields c(v) = 0(1), but then the 1 in EclJ4.3) 

hjr?l'^^ and then we obtain, from Eq.(3.34|, 
4q!i and h = — 4\/2(l + i)a\. Therefore, we 



cannot be cancelled either. Finally, let us try with lixv ^ nn + a/n 
that c(i/) = 0(^1/2) = 0(i^i/2) and from Eq.flol it follows that: a 



have that the highly-damped electromagnetic QNM frequencies are given by 



m 

y 



2n 



23/2^3/2 



o 



(4.4) 



This form agrees with the form of the numerics of Eq.(4.2) (considering that only the leading order in the imaginary 



part is shown) and Fig 12 shows that it also agrees with the numerical data in 



V. HIGH-FREQUENCY BC RESPONSE TO A PERTURBATION 



In this section we consider a perturbation to a Schwarzschild black hole and we investigate the high-frequency BC 
contribution of the £-mode to the response to such perturbation. If the perturbation is given by some initial conditions 
u^i^r'^) = ueXr'.^,t' = 0) and u^i{r'^) = dtUe{r'^,,t' = 0), then the full £-mode response is given by 

/oo 

dr: [Gr*(r,r';t)u-(r:)-H<K)5tGr*(r,r';t)] (5.1) 
-OO 

while the BC contribution to the response is given by 

uf''{n,t)^ / dr:[Gf'^{r,r';t)uT{r:)+uT{r:)dtGf'^{r,r';t)] (5.2) 
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X^' >■• ■••"• :V 



WjJp • , " • ^ 



FIG. 12: Comparison of the QNM frequencies iui„ given by the asymptotic expression Eq.( |4.4[ ), which we will denote by uji„, 
and the numerical data in [35] (after complex-conjugating it in order to account for the different definition of the frequency), 
which we will denote by ijjI^„. For s — 1 and £ — 1. Figs.(a)-(d) are log-plots. Horizontal axis: overtone index n. (a) In blue: 
|Re(tj;'^)|; in green: |Re(a;f„)|. (b) In (dark) blue: |lm(d;[5i) -I- ^j; in green: the leading order for |lm(ii;f„) -I- ^j, i.e., in 
light-blue: |lm(a}"„) + f |. (c): 'Relative error' | |Re(aJiJj)/Re(ti;f„)| — 1|. (d): Same as (c) but for the imaginary part, (e) Plot 
of ■ Re (aj;*^ — a;°„ ). (f) Plot of n^^^ ■ lm(Lu"„ — a;f„). Figs.(e) and (f) show, respectively, that the difference between the 
asymptotic Eq.(4.4l and the numerical data in |39] is only at order for the real part and n~^^^ for the imaginary part. 



We will consider the case u^['{r'^) = dr'^uf{rl), corresponding to an initial wave-packet moving towards the left. 
In addition, for convenience, we take a normalization given by J^^x,'^^* \^Ti''''*)\ ~ ^- define Fe{r,uj) = 



fi{r,ui)e "^''•^ and by analogy with define Fi{r,—iv) — — sm{2T:v)F{r, —iv). Introducing Eq.(2.9l into Eq.(5.2) 
with these initial conditions, we obtain: 



uf^{n,At)=l dv I{r,v), /(r,i/) = -e- 



(5.3) 



+CSO 
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We have defined the function Ct{v) so that, when evaluated at the quasinormal mode frequencies, it corresponds to 
the 'excitation coefficients' Cin defined below for the QNM series. 



We will now investigate the behaviour for large-i^. From Eq.(3.30), 



Fi{r, —iv) ^ c{v) + e 



sin(27rp) 
2-\/p sin(27ri') 



-2vr, 



0,2 

: 1 



(5.4) 



for V 



oo. From Eqs.(3.27) and (3.361, we have that as — ?► oo 

(-l)'*/22cos(27ri7) 



0,2 



1 



\W±\' 
-2vq{v) 



iy[l + 3cos2(27ri^)]' 



(5.5) 



j/3/2 31^(2711^) 



1 



21/ siii(27ri^) 

In order to study the large-i/ asymptotics for Ct{v), we first re-express it as 

/•oo 

Ci{v) = -2vc{v) / dr, e'"-*<(r,) 



(5.6) 



after integration by parts, where we have used Eq.(5.4) and we have assumed that the initial conditions are such that 
the boundary terms are zero, i.e., that F({r^,)e'^^'^lL^Jr^,) — >■ as — > ±oo. 

It is clear that if the initial conditions have co mpa ct support, with u^'^^r^) = for all > i?*, then_|C^(z^)| < 
2 |c(i^)| j^e"^* . Together with the leading-order Eq.( 5.5 1, this determines that the integrand /(r, i^) inEq.(5.3) goes, at 
most, like e'^(~'^*+l'''l+^*) as — )■ oo (ignoring powers of D). Therefore, if the initial data is of compact support, the 
i/-integral in the BC contribution to the perturbation response will converge after a certain time: for At > \r^ \ + R^,. 

Let us now consider the case of non-compact initial conditions, specifically the case of a Gaussian distribution 
centered at = xq and moving towards the black hole, which has frequently been used in the literature: 



1 



/27rcr 



exp 



{r* - a^o) 
2^2 



(f * - Xq) 



(5.7) 



for some cr e M. Note that Eq.(5.7) is the case of Fig. 2 of Leaver [2^ with a = 2 ^^^{2M)^ (apart from a different 
overall constant factor). The high-frequency asymptotics of C^(j^) for the initial conditions (|5.7[) are then given by: 



s = 0,2 

s = 1 



rh 



oo. 



(5.8) 



rh 



In Fig 13 we plot |C£(j^)| as a function of Mv. The leading-orders in Eqs.(5.5l and (5.8) together yield that the 



/2 iy{At |r,| xq) as -> OO (ignoring powers of i/). This asymptotic behaviour 
implies that the z^-integral in Eq.(|5.3) will not converge for any given values of At and r*. See Fig 14 for a plot 



integrand I{r, v) is of the order of e" 



of \l{r,v)e''^\. Leaver's Sec. III. A \2j considers the particular initial perturbation Eq.(5.7) that we have used here. 
Eqs. 54-56 [2], however, only investigate the 'late-time response', i.e., they are obtained via Eqs. 38-44 |2;, which are 
obtained for f ^ 1 and v -^X. This is probably the reason why the large- divergence observed here went unnoticed 
in [2]. This large-j> divergence for any At and r* for Gaussian initial data is in contrast with the convergence for 
At > \r^ \ + in the case of initial data with compact support that we have seen above. 

We expect a similar behaviour for the highly-damped QNM's in the overtone n-sum for the QNM contribution to 
the perturbation response: 



u 



QNM 



{r,,t) = 2j2 



QNM 



{r,,t), 



u 



QNM 



{r.,,t) EE Re 



n=0 



,C,„^>(n,w;„)e-^"'"(^*-'--' 



(5.9) 
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where Cin 



= Ceii'in), the QNM 'excitation factors' are defined by B, 



/{(jJinOHn), and ain is defined via 



as u) 



The n-sum in Eq.(5.9l is over all QNMs in the fourth quadrant of the complex-w 



plane. We note that the radial function gi does not appear in Eq.(5.9) because we may replace it by fi/A°"* at a 
QNM frequency: the two quantities are equal when u — ujin, as follows from the fact that A"^^ 
boundary conditions (2.3) and (2.4). 



and from the 



Comparing Eq.(5.9) with the BC contribution Eq.(5.3), we can say that: iuJin, 2ReX] 



Cpn in the 



QNM contribution 'play the role' of, respectively: dv, ^^^^^^^i , Ci{h') in the BC contribution. This is particularly 

true in the high-damping limit, n oo. 



Let us find the large-n asymptotics of B^n- By comparing Eqs.(2.4) and (3.30) for r — >■ oo we have that A 



s = 0, 2 : A? 



ution for > 0. From Eq.(3.34) it then 



(5.10) 



s = 1 



A", 



as — )■ oo, which is valid since A°^'^^ is the coefficient of the dominant so 
follows that, to leading order as — >■ cxd, 

2-yPsin(27rP) ' 

in the 4th quadrant of the complex-o; plane. In order to compare with |12j . let us note that our definition of tortoise 
coordinate r* equals that in Eq.B.l [12 plus the constant 'm\ Therefore, our coefficients A^^^ and A"^^^ correspond, 
respectively, to l/T^ and e^'^^'^' ^ Rn / T^ , where Rn and are the reflection and transmission coefficients defined 
by Neitzke via Eq.2.10 [T2j. Eq.(5.10) for s = and 2 then agrees with Eqs.2.17 and 2.18 [H], after also taking into 
account a different sign in the dcflnition of lo and converting expression (5.10) from the 4th quadrant into the 3rd 
quadrant by using the sym metrie s (2.5). Note that the different constant of integration '«7r' in = r*(r) is probably 

and 2 differs by a factor e"^'^*" from Eq.A.2 [50]. In the Appendix [a| we 



the reason why A"^^^ in Eq.(5.10) for s 

calculate the relationship between the coefficients A°^^ and A™^ following a different contour in the complex-w plane, 
thus providing a check of Eq.(5.10). We can now calculate that as n — > oo 

s = 0, 2: a^„--37rr,i, B, 
2irfin 



In 



S = 1 



A2 



Bin 



(_l)«+./2j 

27rn3/2 ' 



(5.11) 



where we have used Eq.(3.36) to obtain agn and Eqs.(4.1) and (4.4) for the highly-damped QNM frequencies. The 
behaviour Bgn = 0{n~^) for s = seems to roughly agree with Fig. 2 [B] (although those results are not really meant 
to be valid for large-n). To the best of our knowledge, the expressions in ( |5.10 ) and (5.11 ) for the s = 1 case are given 
here for the first time in the literature. 



The asymptotics as rt — ^ oo for i^£(r*,a;;„) and C^n are those of Fi(r^,—iv) and Ci(y ) in Eqs.(5.4) and (5.6), 
respectively, with the replacement v — > ia;^„. Combining the asymptotics of Eqs.(5.4), ( 5.10| ) and (5.11 ) we find that, 
for n — oo. 



0,2 



1 



u 
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(-1)"A3V2^ 
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iCin (c{iyin)t 
'^tn {c{vin)e~ 



-iy,„(At-r.) _|_ g-i/i„(At+r.) 



(5.12) 



■ui„{At-r,) 



-iy,„(At+r.) 



where a right-moving and a left-moving wave in the radial direction can be seen, both exponentially-damped with time. 
In particular, if the initial data is of compact support vanishing for > i?*, then, for large-n, |Cf„| < |c(;^i„)| ne"'^*^'^ 
and it then follows from Eq.(5.12) that the QNM n-sum in Eq.(5.9) will be convergent for At > -f R^,, just like 
the corresponding BC contribution. 



Let us now look at the case of Gaussian initial data, Eq.(5.7); from Eq.(5. 

(-l)^/2+i+"y3n 



we find as n — > oo 



0, 2 : C£n 

Cin 



S = 1 



2rh 
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(5.13) 



27rAr/i 
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The divergence of C^n for large-n is not cancelled out by any other quantity in Eq.(5.12|, so that the QNM n-sum for 



the non-compact initial data (5.7 



Specifically, from Eqs.(5.12) and (5.13), 



is not convergent, just like the z^-integral in the corresponding BC contribution. 



s = 0,2 
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QNM 



(5.14) 



)/2. 
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In 3 /"a^N 
An 



rhn 



-At + n+x-o)] + e"'^^^*+'*^/'sir 
(-l)"yri^_„ 



An 



At — r^, + xq 



/2n\ 



g ,.(At-f.)/2 g-n(At+f.)/2 



A naive attempt at the calculation of the corresponding large-frequency divergence in the BC contribution indicates 
that it does not cancel out the large-n divergence in the QNM contribution in the case of Gaussian initial data, 
although it is hard to be definitive given the numerous simultaneous asymptotic limits and integrals involved. On the 
other hand, the full perturbation response Eq.(5.1 1 is known to be well defined, therefore, we expect the high-frequency 



arc contribution together with the BC and QNM contributions to be regular, with the divergences in the different 
contributions cancelling each other out. We also expect a similar cancellation between the high-frequency divergences 
from the BC, QNM and high-frequency arc contributions in the case of compact initial data for At < \r^, \ + xq. 




FIG. 13: Log-plot of jC'£(!^)j in Eq.(5.3l as a function of uM for s = Q and 1=1. The blue curve is obtained with the method 

in [29]. The overlapping curve in green is obtained using the large-!/ asymptotics of Eq.( 3.30 1 for in Eq. (5.31. The coefficient 
C(_(v), \lv, for both curves has been obtained by integrating using the built-in function NIntegrate from r = rn up to oo in the 
computational software program Mathematica. 




FIG. 14: Log-plot as a function of vM of \I{r = lQM,v)\ in Eq.(5.3l but without including the factor e "'^ . For s = and 
1 = 1. The darker blue curve is obtained with the method in [29] and the light-b lue cu rve is the inte rpola tion of this data. The 
(partly overlapping) green curve is obtained with the large-;/ asymptotics: Eq.(3.27l for qiy), Eq.(3.30l for Fi, Eq.(3.36l for 



the Wronskian and Ci{u) as in the green curve in Fig 13 
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VI. CONCLUSIONS 



The branch cut contribution to the Green function in the Schwarzschild spacetime has scarcely been investigated in 
the Uterature, except in the small-frequency regime. The challenging technical difficulties for the calculation away from 
that asymptotic regime is probably a reason for it. However, in order to have an understanding of the full response 
of the black hole to an initial perturbation or, in general, to know the Green function globally, we also require the 
non-small- frequency contribution. In [41], the self- force was calculated in the particular case of a static particle in 
(a static patch of) the Nariai spacetime. In such spacetime the 'retarded' Green function possesses no BC (since the 
radial potential falls off exponentially near the 'cosmological horizons'); however, in the Schwarzschild spacetime the 
Green function does possess a BC, the knowledge of which is required if we wish to perform a similar calculation of 
the self-force to the one performed in [41^ in Narai. In this paper we have derived the large-frequency asymptotics for 
the BC contribution to the Green function modes for spins 0, 1 and 2 in the Schwarzschild spacetime - see Eq.( |3.39| 
for the leading order. These asymptotics lead to a divergence of the BC Green function at the 'very early' times 
At < \r^. \ + |r^| and they lead to a convergence of the z/-integral in the upper limit when At > -I- \r'J. In Secjv] 
we have investigated the contributions of the high-frequency BC and the high-frequency QNM to the perturbation 
response given an initial perturbation. We have shown that if the initial data is of compact support within r.^, < i?*, 
then both contributions separately converge when At > |r«| -I- i?*; however, when the initial data is a non-compact 
Gaussian distribution, both contributions diverge separately for any fixed time At and radius . We expect that the 
divergences from all contributions to the Green function (i.e., from the QNM series, the BC and the high-frequency 
arc) to cancel each other out so that the full perturbation response is finite - we leave this study for future work. 



We have also obtained in Eq.(4.4) the highly-damped QNM frequencies for spin-1 (and reproduced existing results 



for the corresponding spin-0 and spin-2 frequencies) in Schwarzschild for the first time in the literature (other than the 
previously- known leading-order for the imaginary part). In order to obtain the leading order of the spin-1 asymptotics 
we had to go up to two orders higher for large-i/ than is necessary for s = and 2. The real part of these spin-1 
frequencies approaches the NIA, unlike for s = and 2, and it does so faster (like n^'^/^) than for s = 1/2 and 5/2 
(which go like n~^/^). 
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Appendix A: Check on the relationship between A'l'^^ and A} 



In this appendix we calculate the relationship between the coefficients A^^^ and A^"^ by following closely the method 
in [12]: we will analytically continue the radial function fi starting at point A in Figjl] then down to point B, where 
arg(r) = 37r/4, then around an argument of '— 37r/2' to reach the anti-Stokes line at arg(r) = — 37r/4, then down 
that anti-Stokes line to radial infinity (i.e., to the 'refiection point' of point A), and finally all around radial infinity 
anticlockwise back to point A (this is Fig.l [H]). We will th en im pose the e xact monodromy Eq.(3.29). This is a 
check on the relationship between A'^^Jj 
do it here only for s — 2. 

We start with the boundary condition /^(r, — zi') 
Eq.(|2!4|). At arg(r) = 37r/4 we have 



and found in Eqs.(3.36) and (5.10) following a different contour. Wc will 



-iv) + AYl^gair, -fw) at the point A in Fig|Tj from 



4 
4 



=V2°^(t) - -3ie*'/2 + e-*'/2 _ _3i5^(r-,z^)e"''-|-ga(r,~ii/)e-''^^ 



(Al) 



The first step follows from Eqs.(3.17) and (3.211; in the second step we have used Eq.(3.25). We solve for gaW^rkiv) 
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and replace in the above boundary condition to obtain 



(0)/ 



(A2) 



We have obtained this expression for arg(r) ~ 37r/4, but it is vahd Varg(r) by analytic continuation. We now explicitly 
evaluate it at arg(r) = — 37r/4. From Eqs.( |3.14 1 together with Eq.(Al) we obtain 



(0), 



4 

4 



(A3) 



at arg(r) — — 37r/4, where in the last step we have again used the top equation in (3.25 1. Introducing these expressions 
into Eq. ( A2 ) we obtain 



/Kr, -iv) ^ ga(r, -iv)e-^^^ [-3^rje"'' + 2^A™ ''J + g„(r, iv)e^^^ [lA^e'^^ - 3A™ ''J (A4) 

at the 'reflection point' of point A in Fig. [l] We can analytically continue this expression all around infinity an- 
ticlockwise and back to point A. We th en equate this expr essio n to the original boundary condition given above 



after applying the monodromy Eq. (3.291, i.e., we equate Eq.(A4) to fi ((r — rh)e^'^^ ,~w) 



^ — liriv A out 



^ AY^^ga{r, Equating the coefficient of ga{r,—iv) gives precisely the relationship between A^^^ and A°^ 



that follows from Eqs.(3.36) and (5.10). Note that the coefficient of ga{r,w) is not to be trusted since when closing 



the contour at infinity it is Re(r*) > and so ga{r,iv) is the subdominant solution there. 
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